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The strong topological insulator in 3D is expected to realize a quantized magneto-electric response, the so- 
called axion response. However, many of the materials predicted to be topological insulators have turned out to 
be metallic, with bulk Fermi surfaces. Following the result of Bergman et al. (Phys. Rev. B 82, 195417 (2010) 
) that the helical surface states of the topological insulator persist even when the band structure gap is closed, 
we explore the fate of the magneto-electric response in such systems. We find a non-quantized magneto-electric 
coupling remains once a bulk Fermi surface opens - a non-universal axion response. More generally we find that 
higher dimensional analogs of the intrinsic anomalous Hall effect appear for every Chern form - non-quantized 
response coefficients for gapless systems, as opposed to quantized transport coefficients in gapped systems, both 
with a topological origin. In particular, the non-quantized magneto-electric response in 3D descends from the 
intrinsic anomalous Hall effect analog in 4D. 

PACS numbers: 



Axion electrodynamics, long confined to the realm of high 
energy physics'! has been recognized recently as a possi- 
ble effective phenomena in solids^. The strong topological 
insulatoi^El at its boundary, is predicted to exhibit the physics 
of axion electrodynamics. The physical consequences of this 
fact are fantastic - ranging from a quantization of magneto 
electric coupling in materials with neither inversion symmetry 
breaking, nor time-reversal breaking, all the way to effective 
magnetic monopole image charges appearing when an exter- 
nal real charge is placed outside the material. 

Perhaps the greatest current issue with 3D topological insu- 
lator candidate materials, which has so far prevented putting to 
the test the interesting predictions mentioned above, is the fact 
that they have for the most part turned out to be metalUcPCni 
The metallic states in these materials do show one character- 
istic feature of the topological insulator - the helical surface 
states. Theory tells us that topological states are only well de- 
fined when the spectrum of the bulk is gapped, yet in these 
materials the surface states coexist with a bulk Fermi surface. 

Having noticed the coexistence of the helical surface states 
and bulk Fermi surfaces, and the observed experimental fact 
that they seem at least mildly robust to disorder effects^ 
Bergman et al. ' ' recently explored the physics of the surface 
states of a hybrid system of metallic and topological insulator 
bands, concluding that it is entirely possible to have helical 
surface states (a litmus test of the topological insulator) coex- 
isting with a bulk Fermi surface. 

The simplest way to think of these "helical metals" (bulk 
metals with helical surface states) is as follows. The minimal 
feature in the bulk band structure of a solid giving topological 
insulator behavior is a point in the first Brillouin zone around 
which the dispersion is approximately that of a massive Dirac 
point - dispersion E = ±-\/ + q^. If a Fermi surface 
opens around some other point in the Brillouin zone, an ap- 
proximate low energy model for this band structure would 
consist of multiple valleys - separate valleys for the massive 
Dirac point, and for the Fermi surface(s) (see Fig. |l(b)| i. Much 
like the case of Graphene (see Fig. 1(a) i, where the two Dirac 
points are to a good approximation independent because of 



the large momentum separation between thenP, the massive 
Dirac point and the Fermi surface(s) will be to a good ap- 
proximation independent when they have a large separation in 
the Brillouin zone. From this simplistic picture we can un- 
derstand that the helical surface states occur at the surface 
as a consequence of the massive Dirac point, and the Fermi 
surface simply does not always interfere with these surface 
states. In particular, the surface states are the Jackiw-Rebbi 
bound states of the Dirac valley. The band structure of Sb, 
expected to be a helical metal" as it exhibits helical surface 
states, can be described in low energies as a number of mas- 
sive Dirac points, in addition to quadratic-dispersing Fermi 
surfaces - this would suggest that the simple multi-valley pic- 
ture for the helical metal applies to Sb directly. 

Using the simplistic multi-valley picture, where the Fermi 
surface and the massive Dirac point are for the most part in- 
dependent, one would naively expect that the axion response 
of the 3D topological insulator would also appear in this (hy- 
brid) system. In this article we explore this question in detail. 
We will find that the axion response can indeed appear, but its 
coefficient is no longer quantized. 

A less obvious, yet even simpler model of the helical metal 
is the Dirac-like Hamiltonian 



-ff = -/-i + gi7i + ?272 + 9373 + Mj5 , 



(1) 



where 71. ..5 are the 4x4 7-matrices, and the chemical po- 
tential is set away from the gap in the spectrum > 
In our convention the 71... 4 are odd under both time reversal 
and inversion, and 75 is even under both, so that this Hamil- 
tonian is time-reversal and inversion symmetric. This model 
qualitatively reproduces most features of the low-energy spec- 
trum of Bi2Se3^ - namely a single massive Dirac point in the 
low-energy bulk spectrum, and a single massless Dirac sur- 
face Fermi surface. Both ([T} and minimal multi-valley mod- 
els can be captured in a generic 4-band model of the form 
H — do(q) + rfa(q) • 7a- Using the 4-band model we will 
show axion response can exist in many gapless systems. 

The axion response action ^axion = j^rf^ J-^ t ■ B = 
8$hc /x t 'l^^i^'^\pdtJ,KdxAp, (from this point on we will set 
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(a) Valleys in Graphene 



(b)Valleys in the helical metal 



FIG. 1: Multi-valley low-energy description of Graphene and of the 
helical metal. Panel [T(b)| shows the two Dirac cones of Graphene, 
and panel 1 1(a)] shows a (dark-shaded) massive Dirac cone at the BZ 
center, and a (light-shaded) quadratic-dispersion Fermi pocket at the 
BZ comer. 



c = h = e = 1) requires the presence of an additional 
field (j> perturbing the system. For systems that are intrinsi- 
cally time-reversal and inversion symmetric, the action must 
remain invariant under these symmetries. The additional field 
(j> therefore can either be odd under time-reversal and inver- 
sion, or if it is a compact variable, e.g. </) = (/) + 20o so that 
it is defined modulo 20o> it may assume the values (p = or 
(j> = (j>Q = —(j)Q. Realizing a spatially varying non-compact (j>, 
that is odd under both time-reversal and inversion symmetries, 
could be accomplished by magnetic impurities, which break 
both, and therefore engender a component that is odd under 
both. With an impurity concentration varying on macroscopic 
length scales, we can realize a slowly varying (j). 

The simplest way to see why we need the auxil- 
iary field (j) is as follows. Naively, one could calcu- 
late a E ■ B response for any given material using lin- 
ear response, which in any system with space-time trans- 
lational invariance would yield S — ^, K^p{x — 

x')A,{x)Ap{x') = /^^^,^^i^,p(g)e^(---')-«A,(x)^p(a;'), 
where x, x' denote space-time coordinates, and q denotes 
a momentum-frequency vector. One needs to expand 
in small q (a gradient expansion) K^p{q) « const + 
(linear in q) + ^fiuXpQf^Qx in order to recover the two deriva- 
tives in dfj_A^d\Ap. For the E ■ B term one would need 



. Unfortunately, q^qx is symmetric 
in its two indices, and such a term would vanish. Therefore, 
linear response cannot give a nonvanishing E ■ B response. 
Remarkably, this also implies that even in a system with inver- 
sion and time reversal symmetries explicitly broken, where the 
term E ■ B does not violate any symmetries, no axion term can 
come from linear response in a translationally invariant sys- 
tem. Consequently, we need the system to be inhomogeneous 
- and this is provided by the auxiliary field 0, varying in space. 
Response theory is easier to calculate assuming a translation- 
ally invariant system, so one can treat the <j) field as an addi- 
tional external field perturbing the system and calculate non- 
linear response instead of linear response. In analogy with 
a gauge field, whose derivatives are the physical observable, 
only the derivatives of (j) appear in the modified Maxwell's 
equations and have measurable physical consequences. 



In the helical metal, the gapless version of the topological 
insulator, we will find that the coefficient in front of the axion 
action is modified to a non-universal value 



, 1 
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(f>E- B , 



(2) 



with c 7^ 1 and depending on the system details. The non- 
quantized magneto-electric response coefficient is analogous 
to the Hall conductivity in the intrinsic anomalous Hall effect 
(AHE) in 2D Ferromagnetic metals - because there is no gap, 
the Hall conductance is not quantized, even though it does 
have a topological origin - the Berry phase - and when opening 
up a bulk gap, the Hall conductance is quantized. 

We will start by demonstrating that an analog of the AHE 
occurs in gapless non-interacting systems in any odd space- 
time dimension d = 2n + 1, and from it, we will understand 
how the axion response appears in the helical metal. Being 
the first attempt at calculating this response, we will ignore 
the effects of disorder and of finite temperature. 

The Chern Simons forms in electrodynamic response in 
d = 2n+ 1 space-time dimensions, ar&Sl 

Schem — Xn^fii...p2n+l / ^Ml "-^^2 • ' • ^^2,1 ^A'2n + l • (-^) 

J X 

The coefficient x„ can be found using response theory 



i{n + I) {2n + 1)1 



2n+l 

n iGdp.G-') 



(4) 



where G = [iutn — H{q)] ^ is the single particle Green's 
function, and du — with q^ the real frequency, and qi 2n 
being the Brillouin zone coordinates. 

We will calculate this coefficient for a Hamiltonian of the 
form H ~ rfo(q) + dai^i) ■ 7a where the matrices 7i...2n+i 
satisfy a Clifford algebra {7a, 7b} = 2Sab, as well as the iden- 
tity 



T'f [7ai7a2 ■ ■ ■ 7a2„7a2 



2meai. 



•0271 + 1 • 



(5) 



With do — this is precisely the form of a minimal model 
for the Chern insulator in d = 2 + 1 and for the time- 
reversal invariant 4D Quantum Hall state in rf = 4 + lISISl 
(which realizes the strong topological insulator when reduc- 
ing to d — 3 + 1). With do closing the gap in some re- 
gions of the Brillouin Zone, this describes the helical metal. 
First, we "normalize" the Hamiltonian by i? = \d\, by writ- 
ing H = Rh = R{ro + raja), where now — 1 and 
g ^ GR ^ R[iujn-Hy^ = [iujn/R~hy^. It follows 
gd^g^^ + dp ln(i?), in which the last term is 



that GdpG-^ 



a scalar, and so can be moved outside of the trace in Q. As a 
result, when expanding the formula (|4]i, any term where more 
than one instance of GdpG~^ — > dp ln(i?) appears, will van- 
ish, since dp \n{R)dp> \n{R) is symmetric in the indices /i, /i' 
and we have the completely anti-symmetric tensor multiply- 
ing the entire term. From the remaining cases where zero or 
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one instances of Gd^G ^ 



ln(i?) appear, we have 



Ml---M2n + 1 



^" i(n + l)(2n + l)! ' ^ 



2n+l 



2n+l 



+ ^ 5,,^ HR)Tr 
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We will show the second term vanishes, by integrating by 
parts with respect to some g^^. We first rewrite gd^^g^^g = 
—d^.g and then integrate by parts. Acting on the scalar term 
outside the trace, we get d^j^^d^j ln(i?), which is symmetric 
in the two indices and thus vanishes at the hands of the anti- 
symmetric tensor Similarly, every instance in which the 
derivative acts on any existing derivative term df^^g^^ also 
yields an expression symmetric in its indices, and vanishes. 
Otherwise, it must act on an instance of g, giving d^^^g = 
—gdfj^.g^^g and producing pairs of expressions of the type 
Tr [. ^..gd^^g~^gd^,g-^ . . ]+Tr [. . . gd^.g^^ gd^^g~^ . . .] 
symmetric in fii,i, and thus also vanishing. There are 2n in- 
stances of g in the trace, of which 2 get absorbed into the 
initial d^-^g, leaving 2n — 2 instances of g on which 9^. acts 
on after integrating by parts. Since this is an even number, all 
such instances pair up as described above, and the whole sum 
vanishes. We are left therefore, with the same expression as 
in (|4]) with g replacing G. This will prove easier to calculate. 
It is worth mentioning in passing that for the case of an in- 
sulator, this is precisely the same procedure of deforming the 
bands to flat bands, used in Ref. 2 

Using the Clifford algebra, the Green's function can be re- 
written as g = ^i'^nZl'^'^-T ^ = 0). The terms d^g-^ 
will have a scalar part d^{iuj/R — tq), which can leave the 
trace. If two or more instances of this scalar part appear, it 
will be symmetric in the two derivative indices, and vanish 
under the anti-symmetrization. Therefore at most one such 
scalar part can appear, but this is also the minimal number, 
because one of the derivatives has to be with respect to fre- 
quency dog = l/R. As a result, for all the other terms, 
we can ignore the scalar part of the current terms and take 
gd^.g^'^ = ^'('^i^fl-°o|2''°i° i-lhbd^n. For every permuta- 
tion of /ii...2n+i we can cyclically rotate the dog term to the 
left end in the trace formula, yielding 



n+l. 



Xn 
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i(n+ l)(2n)! 



gdog ^ Yl (99^9 ^) 
3=1 
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where now the indices run only over momentum coordinates. 
From = 1 we get f ■ dfj^f = 0, which we can use to find 

gd^g-^gd^g-^ = [^'rJ/tlf)!^-"!] - Plugging this in Q, we 
arrive at 



Xr, 



(-1) 



n+l. 



i{n + l){2ny. X 



[(zw/i?-ro)2-l]"i? 
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We first turn our attention to the trace over 7 matrices. Due to 
the anti-symmetry in the indices, the product of jai 2n 



should also be anti-symmetric in its indices. In particu- 
lar, those indices should be all distinct. It is then easy 
to show that the trace over these 2n matrices must vanish, 
Tr [7a J . . . 7a, ,J = 0. Together with (|5]l we then find that 



(-1) 



n+l. 



(n + l)(2n- 1)! 
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Performing the T 



Matsubara frequency integral, we get 

_ (-i)"+'7^r(«+l/2) ^., Ir..n-.vhcrc 

-00 [(,<^_,.o)2_i]"+i - r(i+«) ^y'^ |ro|jwnere 

r{n) is the Gamma-function, is a step function, and we 

absorbed R > into ui. Using this, we find 



<^ai...a2„+i 



2n 



2n 

•^9 j = l 



(10) 

where (?i...2n run over the Brillouin zone coordinates. Only 
those areas of the Brillouin zone (BZ) where there is a lo- 
cal gap Id I > I do I contribute to the coefficient %„. If there 
is a gap everywhere in the BZ, the integral gives a quantized 
value. The quantization can be understood by realizing that 
the integrand calculates the (directed) infinitesimal hypersur- 
face element, and the integral measures the total hypersurface 
swept by the mapping r(q) from the BZ to the unit hyper- 
sphere S2n (the 2n-sphere). This integral will give an integer 
times the hypersurface of the 2n-sphere ^ 



ing us the transport coefficient Xn — ^ (n+i 



2n 



(2n)! 



glV- 

7 — rTYm-viT, where 
N is the topological integer For n = 1 this describes the 
Hall conductivity, which in the gapped case gives the integer 
quantum Hall effect, and in the gapless case gives the intrinsic 
anomalous Hall effect contribution'^. For n = 2 the gapped 
case gives the 4D quantum Hall effecfSIS! while for the gap- 
less case we get a non-universal transport coefficient. The 
n — 2 Chern response is the precursor of the axion response 
in 3D, and so our results strongly suggest that such a response 
can be found in 3D gapless systems, to which we turn now. 

The simplest way to think of the axion response in 3 di- 
mensions is to start with the the n = 2 Chern form ([3]l, 
and replace the extra vector potential component A4 with 
the auxiliary field </), mentioned earlier, which couples to 
matter fields via the free-fermion propagator G = G{q, (j)). 
Note that since A4 is odd under time reversal and inver- 
sion, (p must be odd as well, fitting its non-compact vari- 
ety. The auxiliary field ip couples to an extra current operator 
J4 = —d^G^^{q, (t>)\tp^Q, and all the other current operators 
are now — — 9^G'~^(q, 0)|0=o- In addition, the 4th spa- 
tial dimension gets dropped, and the anti-symmetric symbols 
now include a 5th coordinate representing the extra field (f). 
The derivation of the non-linear response follows through in a 
similar manner as described above for the d = 2?i + 1 case, us- 
ing G{q, (f>), except perhaps the integration by parts arguments 
used when converting to the "normalized" bands. However, 
the integration by parts used on (|6]l, could be done with re- 
spect to any one of the 2n > 2 different spacetime indices /i^, 
so we can easily choose one that is a frequency -momentum 
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FIG. 2: Calculation of the magneto electric response coefficient in 
a model system l[T2j realizing both a gapped and gapless spectrum. 
The transport coefficient c is quantized as a long as a gap remains 
l/^l < M. Once the gap is closed, c is nonquantized. 



component rather than cj). The calculation therefore follows 
through exactly in the same manner as in d = 4 + 1, and we 
get after a bit more algebra the coefficient from (|2]) 

c = 67r^eai...a5 / Q{\d\ ~ IdoDraidq^ra^dq^ra^dq^ra^d^ra^ 

Jq 

(11) 

To demonstrate that c can be nonzero in a gapless system, 
we will calculate it explicitly for ([T]|, coupled to (j> as follows 

H{(f>) = -/i + gi7i + 9272 + 9373 + (/'A/74 + Af75 . (12) 
Using the formula ( fTT) , we can easily find 

SAPq^dq 



(g2 + Af 2)5/2 



e{q' + hp - M^) 



(13) 



Calculating numerically for a range of values for the chem- 
ical potential ^ we find the results of Fig. [2] which clearly 
indicate a quantized value c = 1 for the gapped case and a 
non-quantized value for the gapless case. 



In summary, we have shown that in higher dimensions an 
analog of the intrinsic anomalous Hall effect exists in gapless 
systems closely related to Chem and topological insulators. 
As a consequence, the 3D descendant of the 4D effect is a 
non-quantized axion response in 3D gapless systems. 

The future directions following what we have shown here 
are interesting and numerous. Borrowing from the what is 
known about the AHE— , we expect that including disorder 
will bring about corrections to the axion response - both from 
a finite lifetime for quasiparticles, and from vertex correc- 
tions. Given that the system is gapless, these coiTections 
may be significant, and one must proceed with caution. Still, 
given the robustness of the AHE in d = 2 + 1 to disorder, 
it would be surprising if its higher dimensional analogs (the 
non-quantized axion response among them) were completely 
washed out by disorder effects. 

Another question that arises is how Coulomb screening will 
affect the electromagnetic system with this strange new re- 
sponse - contrary to the topological insulator, the helical met- 
als will have free carriers in the bulk that can screen electric 
fields. We leave the resolution of these questions to future 
work. The combined issues of disorder effects and screen- 
ing make it difficult at this stage to predict how one would go 
about measuring experimentally the axion response in gapless 
systems. However, probing AC response (where screening is 
no longer effective) may reveal the axion response. We remain 
hopeful that the non-quantized axion response could appear in 
helical metal materials doped with magnetic impurities, per- 
haps in Bi2_2;Mna.Te3'^, doped to form a bulk Fermi surface. 
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